Tunneling mediated by conical waves in a ID lattice 
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The nonlinear propagation of 3D wave-packets in a ID Bragg-induced band-gap system, shows that 
tranverse effects (free space diffraction) affect the interplay of periodicity and nonlinearity, leading 
to the spontaneous formation of fast and slow conical localized waves. Such excitation corresponds 
to enhanced nonlinear transmission (tunneling) in the gap, with peculiar features which differ on 
the two edges of the band-gap, as dictated by the full dispersion relationship of the localized waves. 
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Localized wave-packets (LW) characterized by a con- 
ical envelope, and in particular X waves, have been in- 
voked as a new paradigm that allows to explain several 
nonlinear optical wave phenomena observed in dispersive 
media, such as fsec trapping via frequency doubling and 
filamentation in Kerr-like media [H, H, H S] ? ^^^^ 
spatio-temporal localization in lattices [5|. Unlike soli- 
tons, these LW exist also in the linear limit where they 
are expressed as superpositions of Bessel beams, well 
known in different context of physics 0, 0, @] • The nonlin- 
earity "dresses" these solutions, playing a key role (often 
strongly dynamical [5, ^3^ ) in their formation and morph- 
ing process. 

It has been suggested that LW can also exist as enve- 
lope modulations of linear Bloch modes in periodic media 
such as ID lattices (with free diffraction in the trans- 
verse direction) tell, or 3D generalizations, namely pho- 
tonic crystals |1Q|.Iii|. In this case they play a major role 
around the forbidden band-gaps where, however, strong 
reflection is expected to severely hamper their formation 
in the linear regime. While in principle the nonlinear- 
ity can be envisaged to favour the energy coupling and 
thus the excitation of LW, the dynamics of this process 
is completely unexplored so far. Furthermore, one can 
argue that conical LW can compete with gap soliton bul- 
lets (GB) [13, H [li, which are the multi-dimensional 
analog (usually under more restrictive conditions) of well- 
known gap solitons of ID settings [lH, observed, e.g., as 
slow light pulses in fiber Bragg gratings |l6l, llTI ]. spatial 
envelopes in array of coupled waveguides [l8|, or atom 
clouds in a ultracold Bose condensed gas with standing 
wave optical lattice [19]. 

In this Letter we investigate the dynamic role of LW 
owing to additional free transverse dynamics (diffrac- 
tion) included in a system with stop-band (whose ID 
lon gitu dinal dynamics and gap solitons are well known 
I5I . Til . IitI . [Til . [Tit ) . Employing excitations easily imple- 
mentable in experiments (unidirectional and bell-shaped 
beams), we show numerically that, in lattices of finite 
length, highly transmissive nonlinear behavior in the 
stop-band is mediated by the excitation of conical LW. 
The features of this process can be understood on the ba- 



sis of the linear multi-dimensional dispersion relationship 
of fast and slow LW, which entails a different dynamics on 
the two edges of the stopband. While on the lower edge, 
LW have a primary role to determine the enhanced non- 
linear transmission, on the upper edge LW compete with 
GB, and their excitation stems from components tunnel- 
ing through the gap in transverse wave- number. Thus 
our results suggest that gap solitons are not the sole en- 
tities that rule multi-dimensional nonlinear transmission 
through a stop-band. 

We start from the following continuous model that 
rules (in dimensionless units) the propagation around a 
single Bragg resonance of a period ID system [ol, [H, [T5[ 

(^idt + id, + i Vi^ + + X (l^+l' + 2 I') = 0, 



(^idt - id, + ivi^ + ii+ + X (l^-l' + 2 \u^\ 



(1) 



u- = 0, 



where t and z are normalized time and distance (along 
the ID lattice), and =9^ + 9^ is the transverse La- 
pacian. Though we expect a similar scenario with radi- 
ally symmetric beams in two transverse dimensions, we 
discuss the computationally less demanding case of one 
transverse dimension (V^ = dy). In this case Eqs. ([1]) 
hold for Bragg optical gratins in a slab waveguide, or for 
Bose-Einstein atom condesates with free kinetic spread- 
ing (no trap in y) transversally to a standing wave poten- 
tial sin^ {kg z/ 2). In the former case y are measured in 
units of To = {VgT)-\ Zo = T-\ Yo = (2A:onor)-i/2, re- 
spectively, whereas u± = \//co|^2/|/2r/£^±, E± are coun- 
terpropagating envelopes at Bragg frequency, ko and 77 
are the relative vacuum wavenumber and impedance of 
the host medium, F is the Bragg coupling coefficient, and 
7121 the Kerr nonlinear index. In the latter case, Eqs. ([T]) 
derive from the mean-field Gross-Pitaevskii (GP) equa- 
tion [9] with To = n/F, = h^kg/(2mT), Y^ = h^/mV, 
m being the atomic mass whereas 4F is the peak energy 
of the periodic potential. = a~^^'^u± are mutually 
Bragg scattered components of the atomic wave- function, 
with a = {Anh^ /m)as^ as being the scattering length. In 
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FIG. 1: (Color online) Spatio-temporal dispersion curves 
kz = 00 /v of LW in the frequency plane ky^u. Outside the 
forbidden band (region between the two parabolas), the LW 
Ught-line (t' = 1, blue curve) divides the regions of subluminal 
LW {v < 1, green thin curves) and superluminal LW {v > 1, 
thick red curves). 



the following we restrict the discussion to the focusing 
(or attractive) case i.e. x = sign{n2i) = —sign{as) = 1, 
though similar phenomena can occur also in the defocus- 
ing (repulsive) case. 

The reconstruction of the overall field or wave func- 
tion (apart from inessential common phase shifts) u = 
Uj^ exp {ikgz/2) + U- exp {—ikgz/2), shows the nature of 
the envelopes u± as modulations of Bloch modes. 

Dispersion and localized wave spectra. — The linear 
dispersion [u± ~ ex.p{ikzZ + ikyi/ — iut)] relationship as- 
sociated with Eqs. ([T]) reads as k^ = {uj - ky/2)'^ — 1, 
which entails a lower (LB) and upper branch (UB) sep- 
arated by a forbidden band in uo (or energy). This gap 
of edges a;± = /c^ /2 ± 1 , has constant width and central 
frequency shifting upward (parabolically) with ky (tilted 
waves with ky ^ ^ see an effectively shorter period). 
While the UB u;-^{ky^kz) has always positive curvature, 
it is near the saddle-shaped (i.e., d'^co/dkld'^co/dky < 0) 
LB LO-{ky^ kz) that X-shaped LW has been predicted on 
the basis of a multiscale perturbative approach Q . How- 
ever, at the linear level (x = 0), more general existence 
conditions for non-dispersive LW require simply a linear 
dependence kz = u;/v^ where the arbitrary parameter 
V = dkz/dix)~^ stands for the group velocity of the LW 
travelling undistorted [2,]. Upon substitution in the dis- 
persion relationship we obtain the dispersion curves of 
LW, namely (^)^ = {u — A:^/2)^, shown in the plane 
{ky^uo) (see Fig. [T]) for various values of v below and 
above the luminal curve v = 1, denoted henceforth as 
the LW light-line. The dispersion curves do not involve 
evanescent waves and lie outside the gap (white region 
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FIG. 2: (Color online) Excitation close to LB {uj = -0.99). 
(a) Transmission Vs input fluence E] (b) Output pulse delay 
T Vs E. The delay is calculated as the position of the first 
peak in the transmitted 2/— integrated pulse profile. 



bordered by parabolas). A wave-packet with spectrum 
along one of such curves represents a LW with well de- 
fined velocity [22[, while quasi-stationary superpositions 
of LW \d\ with defined mean velocity result from spectra 
well concentrated in the different regions of Fig. [H In 
the nonlinear regime, as we show below, the nonlinearity 
acts as the triggering mechanism of the spectral reshap- 
ing of the input into such type of LW. We point out that, 
though LW exist above the UB edge where they corre- 
spond to O — waves [7], for excitation in the band-gap 
(regardless whether in proximity of the LB or UB) those 
playing a role are those below the LB edge, which are of 
the X type [7] due to the opposite curvatures. 

Numerical method. — We consider the boundary 
value problem given by Eqs. ([T]) subject to unidirectional 
(forward) gaussian excitation 



ii+(t, X, z = —L) — uq exp 



ex.p{—iujt) 



(2) 

and U-{t^y^z = L) = 0, where co is the frequency 
(energy) detuning that fixes the position in the gap, 
and E = tti^o^o^q input square norm (fluence). 

Though we have performed extensive numerical runs, we 
show below a typical case obtained for 2L = 4, = 30, 
to = 12, yo = 6. Since Eqs. ([1]) and ([2]) cannot be 
solved by standard beam-propagation tecniques [lO] (the 
main difficulty arises from the presence of a finite in- 
tegration region in the propagation direction z and an 
infinite transverse domain in ?/), we adopted a pseudo- 
spectral technique based on a Chebishev-collocation ap- 
proach along 2:, in conjuction with second-order split-step 
(evolution variable t, Fast-Fourier trasform in y) [21]. 

Lower branch excitation. — In Fig. [2^ we show the 
transmission T for an excitation inside the gap in proxim- 
ity of the LB. At low fluence E there is a residual trans- 
mission (owing to the finite size of the system), which 
decreases at moderate fluences consistently with the non- 
linear shift of the gap toward lower frequencies (energies). 
However, at sufficiently high fluences, the trasmission in- 
creases with E. This nonlinear self-transparency cannot 
be ascribed to stable GB, which do not exist close to the 
LB [13] (on the LB gap soliton are unstable even in the 
plane wave limit = 0). However, the phenomenon can 
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FIG. 3: (Color online) (a,b) Output spectrum of at (a) low 
fluence E — 1^ (b) high fluence E — 12{) (inset: zoom around 
ky = 1). The red dashed line denotes the forbidden band, 
the blue thick line is the LW light-line. (c,d) corresponding 
output spatio-temporal profile for (c) ^ = 1; (d) E — 120. 
The vertical dashed line gives the reference instant of the peak 
of \u+\^ at ^ = 1. 
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FIG. 4: (Color online) As in Fig. [2]for = 0.99 (LB). 



shown in Fig. [3]i. 

Upper branch excitation. — The excitation in the 
proximity of the UB is characterized by the existence of 
GB which provide the main tunnelling mechanism [13]. 
Evidence of this process is reported in Fig. [4^, showing 
that the trasmission T increases abruptly around £^ 40. 
Above this threshold, T remains high, while it exhibits 
clear peaks. In addition, as shown in Fig. [Dd, the delay 
versus fluence E displays a cusp in correspondence with 
the threshold. This implies that the generated nonlinear 
waves have a decreasing velocity for increasing fluence up 
to the threshold value, after which this trend is reversed. 



be explained by looking at the spectra and t — y snapshots 
(top view) of the fields at the output {z = L). In Fig. [3] 
we compare the case of low {E = 1) and high {E = 120) 
fluence. While the output ky — uo spectrum at £^ = 1 is 
well concentrated in the gap around small wave-numbers 
(Fig. [3^), at high fluence the strong self- focusing causes 
the spatial spectrum to widen and fill the region of large 
wave- numbers (Fig. [3)3). As also evident in the inset, the 
energy is progressively concentrated in proximity of the 
LW light-line, with a clear peak in the region of fast LW, 
arising from tunneling in wave-number outside the band- 
gap. Correspondingly the t — y snapshots show an almost 
undistorted (nearly gaussian) beam at low fluences (Fig. 
[3b), which at large fluences, develops clear conical tails 
along the trailing edge (Fig. [3li). Importantly the lead- 
ing edge of the pulse arrives earlier in the nonlinear case 
consistently with the spectral peak in the region of fast 
LW. However, also spectral components in the region of 
slow LW are present as clear from the inset in Fig. [3)3. 
The interplay between these two components explains the 
behavior of the output delay vs. E displayed in Fig. [2)3. 
As £^ > 50, the enhanced transparency of the structure 
is accompained by a growing (with E) delay. This cor- 
responds to the excitation of slow LW components (see 
Fig. [3)3). However we also found that energy is coupled 
to frequencies in the region between the lower band-edge 
and the LW light-line (see inset in Fig. [3)3). The latter 
process becomes dominant at large E^ being responsible 
for the delay to decrease with for £^ > 100. In the 
spatio-temporal domain, the fast LW corresponds to a 
bell-shaped leading edge of the transmitted pulse, while 
the slow LW can be associated to slow conical tails, as 



This dynamics is a direct consequence of the multi- 
dimensional nature of the problem and the existence of 
LW solutions. In Fig. [5^1, d we show the ky — uo spec- 
trum and (^, t) snapshot at the output z = L, in corre- 
spondence of the first peak in the transmission T. As 
shown, the tunneling wave-packet has a spectrum which 
remains well confined inside the gap around low wave- 
numbers. Correspondingly the beam, which has under- 
gone self- focusing in the transverse variable, has substan- 
tially bell-shaped section (elliptical top view) with very 
weak V-shaped tails. Such evolution indicates that the 
transmission peak is basically associated with the forma- 
tion of a GB tunneling through the lattice. The situa- 
tion remains nearly unchanged until the second peak of 
transmission at £^ = 64 is reached. Here, the situation 
is radically different as clear from the spectrum in Fig. 
[5)3. Owing to strong self- focusing, a component of the 
spectrum has clearly tunneled (mainly horizontally, i.e. 
in ky) through the gap to the LB. In t^y domain this 
spectral feature has one to one correspondence (as can 
be shown by filtering and anti-transformation) with the 
strong enhancement of the V-shaped tails characteristic 
of X waves (see Fig. [5^). Since the LB spectral compo- 
nent lies in the region of fast LW (i.e., the tiny region 
between the LW light-line and the band-edge), we con- 
clude that the transmission peak is associated with the 
concomitant excitation of a such type of waves. The fast 
character of the LW of the X type is also clear from Fig. 
[5^, which shows that the tails emanate from a wave com- 
ponent that anticipates the slow GB component. The 
tunnelling process due to the generation of the fast X- 
wave goes on until the fluence is increased up to the next 
transmission peak at £^ = 90 (see Fig. [4^), where a sec- 
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FIG. 5: (Color online) Output features of for UB excita- 
tion {uj = 0.99). Output spectra (a) E = 40, (b) E = 64, (c) 
^ = 90 (zoom around LW light-line); Corresponding inten- 
sity (t,y) profiles: (d) E = 40, (e) E = 64, (f) E = 90. The 
dashed line in (e,f,g) gives the reference instant corresponding 
to peak of \u+\^ at ^ = 40. 
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FIG. 6: (Color online) Level plot of spatial intensity \u+\^ + 
\u-\^ for UB excitation [uj — 0.99), at two instants: (a) t — 
40, both the fast conical LW and the GB are trapped in the 
lattice; (b) t = 46, the fast LW has left the lattice and only 
the slow GB remains. 



ond LW is generated, as shown by the splitting of the LB 
spectral component (see Fig. [Sb), and the more struc- 
tured fast wave component in Fig. [5f. The excitation 
of a fast LW explains the cusp in the delay time r in 
Fig. [Dd: as energy is increased the bullet develops fast 
LW tails that eventually precede the GB and correspond- 
ingly (cusp in Fig. [4)3) the delay starts to decrease with 
E. This process is even more evident in Fig. [6l where 
we show the intensity profile for £^ = 50 and 2L = 8 (a 
longer structure makes the process more evident), at two 
different instants. For t = 40, the intensity distribution 
clearly unveils a X-shaped pulse which is in proximity of 
the output. At t = 46 the fast LW has left the structure, 
and only the slow GB remains. 

Conclusions. — In summary the study of transverse 
effects a ID Bragg lattice unveils the role of fast and 
slow LW, characterized by specific spatio-temporal spec- 
tra. Their excitation foster a variety of self-trasparency 



effects, either in conjuction with multi-dimensional gap 
solitons, or as an alternative to them. Measuring the 
spatio-temporal spectra or the trend of the output delay 
with input fluence could be used to monitor the excita- 
tion of LW. Our results have implication in all of the 
field where multi-dimensional solitary wave and period- 
icity play a role, and more specifically in Bose Einstein 
condensation and nonlinear optics. 
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